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ABSTRACT 

Context. Core condensation is a critical step in the star-formation process, but is still poorly characterized observationally. 

Aims. We have studied the 10 pc-long L1495/B213 complex in Taurus to investigate how dense cores have condensed out of the 

lower-density cloud material. 

Methods. We have observed L1495/B213 in C 18 O(l-0), N 2 H + (l-0), and SO(J N =3 2 -2!) with the 14-m FCRAO telescope, and com- 
plemented the data with dust continuum observations using APEX (870 /jm) and IRAM 30m (1200 /im). 

Results. From the N 2 H + emission, we identify 19 dense cores, some starless and some protostellar. They are not distributed uni- 
formly, but tend to cluster with relative separations on the order of 0.25 pc. From the C ls O emission, we identify multiple velocity 
components in the gas. We have characterized them by fitting gaussians to the spectra, and by studying the distribution of the fits in 
position-position-velocity space. In this space, the C 18 components appear as velocity-coherent structures, and we have identified 
them automatically using a dedicated algorithm (FIVe: Friends In Velocity). Using this algorithm, we have identified 35 filamentary 
components with typical lengths of 0.5 pc, sonic internal velocity dispersions, and mass-per-unit-length close to the stability threshold 
of isothermal cylinders at 10 K. Core formation seems to have occurred inside the filamentary components via fragmentation, with a 
small number of fertile components with larger mass-per-unit-length being responsible for most cores in the cloud. At large scales, 
the filamentary components appear grouped into families, which we refer to as bundles. 

Conclusions. Core formation in L1495/B213 has proceeded by hierarchical fragmentation. The cloud fragmented first into several 
pc-scale regions. Each of these regions later fragmented into velocity-coherent filaments of about 0.5 pc in length. Finally, a small 
number of these filaments fragmented quasi-statically and produced the individual dense cores we see today. 
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1. Introduction 

Dense cores are the sites of individual or binary stellar birth 
( |Benson & Myers|1989[|di Francesco et al.|20 07 ), and their con- 
densation from the ambient cloud represents a critical step in the 
process of star formation. Core formation, however, appears as 
a bottleneck in the road from clouds to stars. At any given time, 
only a small fraction of the material in a cloud is in the form 

2007} , and this ap- 



of dense cores (< 10 %, i.e., Enoch et al 



parent difficulty making dense cores is likely connected to the 
low efficiency of the star-formation process (Evans 2008 ). Core 
formation, in addition, seems to play a key role in determining 
the mass of the final star, since the distribution of masses among 
starless cores in a cloud mimics the initial distribution of stellar 
masses ( |Motte et aL|l998[ |Alves et aL|2007) . Core formation, 
therefore, represents a crucial gas transition by which a cloud se- 
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lects a small fraction of its material to form the next generation 
of stars and leaves the rest as a sterile remnant to be dispersed 
into the more diffuse interstellar medium. 

Despite its critical role in star-formation, the process of core 
formation is still poorly understood. A number of condensation 
mechanisms have been proposed over the years, ranging from 



quasi-static contraction mediated by ambipolar-diffusion (Shu 
|et al.||1987] |Mouschovias & Ciolek| [l999) to gravo- turbulent 



fragmentation driven by su personic motions (|Padoan et al. 



2001; |Klessen et aL]|20()5| Vazquez-Semadeni et al. 2005). 



Observations, however, do not favor clearly any single model, 
partly due to the intrinsic difficulty in measuring time scales 
of core evolution and magnetic field intensities, which are the 
main parameters that distinguish the different models (Ward- 
|Thompson et al.|2007 for a review). 

New approaches to study core formation are needed, and a 
promising strategy is to analyze in detail the velocity field of 
the gas, comparing the internal motions of the dense cores with 
those of the surrounding low-density gas from which they have 
condensed. Following this approach, Hacar & Tafalla ( 201 1\ 
found that the gas surrounding the cores of the L1517 cloud 
is subsonic and quiescent, like the gas inside the dense cores, 
in contradiction with the predictions from models of gravo- 
turbulent fragmentation. More interestingly, the lower-density 
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gas in L1517 forms a network of filaments whose velocity fields 
are continuous and subsonic over scales of about 0.5 pc, which 
is significantly larger than the typical core size (« 0.1 pc). These 
so-called velocity-coherent filaments seem therefore to consti- 
tute the parent structures from which the cores form. The ob- 
served continuity between the large-scale velocity field of the 
filaments and the internal velocity gradients of the cores indi- 
cates that the transition between the two regimes involves little 
kinematical changes, and in particular, an absence of shock com- 
pression. 

The analysis of L1517 suggests that core formation in this 
Taurus cloud has occurred in two steps, with the velocity- 
coherent filaments forming first and the cores fragmenting later 
from the already-quiescent filament gas. In this scenario, turbu- 
lent dissipation precedes core formation, and gravitational frag- 
mentation of the filaments is the final step in the core-formation 
sequence. 

Independent work using dust continuum data from the 
Herschel Space Observatory has revealed that filamentary struc- 
tures in clouds are ubiquitous, and that filament-based core for- 
mation is likely to be a widespread process in both low and 
high mass star-forming regions ( Andre et al. 2010; |Molinari| 
et al. 2010| [Ar zoumani an et al.||2011|i . (See also Schneider &| 



Elmegree n"T|1979| |Johnstone & Bally| [Y999} 



Myers||2009| for previous work emphasizing the importance of 



Hartmann 2002 



filaments in core formation.) Understanding core formation in- 
side filaments has therefore become an urgent task, and kine- 
matic information is a needed element to complete the picture 
of filament physics. For this reason, we have selected the most 
prominent filamentary region of the Taurus cloud, which consists 
of the B211, B213, B216, B217, and B2 18 dark pa tches ( |Barnard| 
1927 1 and ends in the L1495 cloud ( Lynds|1962 l, and we have 



subjected it to observations similar to those used to study L15 17 
Our cloud of study, which we will refer to as L1495/B213 for 
brevity, appears in optical images as a several-degree-long fila- 
ment first noticed by [Barnard ( 1907| l, who incidentally used this 
and other Taurus filaments to argue that at least some dark ob- 
jects result from the obscuration by a physical substratum, and 
not from the mere absence of background stars as previously 
thought. Due to its striking appearance, the L1495/B213 com- 
plex has been investigated and characterized by a number of au- 
thors over the years. Initial extinction maps were presented by 
Gaid a et al~| ([1984) and Cernich aro et al.|(|1985|l u sing optical 
data, and more recently by |Schmalzl et al.| ( |2010) > using deep 
NIR observations. The gas component of the cloud has been 
studied systematically in 13 CO and C ls O by |Duvert et al.| ( |l986| , 
who noticed a complex velocity structure that was interpreted as 
resulting from colliding filaments. Additional large-scale maps 
of the region in CO isotopologues have been present ed by|Heyer| 
et~al~| (fT987)l, |Mizuno et al.| ( fT995> , [Onishi et al.| (|1996), and 



Goldsmith et al.| ( |2008| l. Very recently, |Palmeirim et al.| ( |2013| l 



have presented high dynamic range Herschel images of the dust 
continuum emission. 

Embedded in the relatively low-density CO-emitting gas of 
L1495/B213, lies a population of dense cores, which has been 
explored by |Benson & Myers | ( [T989) and |Onishi et al.| ( |2002| l us- 
ing NH3 and H"CO + observations, respectively, with addi tional 
and more limited N2H + observations presented by Tatematsu et 



[aL]([2004). Some of these dense cores are starless, while other are 
associated with young stellar objects (YSOs) of different ages. 
This population of YSOs has been the subject of a number of 
dedicated studies, most recently by Luhman et al. ( 2009 \ and 



Rebull et al.| ( 2010 >, and by the dedicated outflow search from 
Davis et al.i poToT 



The above observations, and additional work reviewed by 
Kenyon et aL] ( |2008| l, show that the L1495/B213 complex has 
been and still is an active site of star formation in Taurus. In 
fact, all stages of the star-formation sequence can be found 
in the cloud, from starless cores that may have recently con- 
densed, to Class and Class I protostars with active outflows 
and signs of accretion, to classical and weak T Tauri stars in 
their path to the main sequence. This rich population of objects 
makes L1495/B213 an ideal laboratory to study core and star- 
formation, with the additional advantage that the region has a 
well-defined filamentary geometry that simplifies its mapping. 

2. Observations 

The main dataset used in this paper consists of observations car- 
ried out with the FCRAO 14m radio telescope during several ses- 
sions between March 2002 and November 2005. In each session, 
the telescope was equipped with the 32-pixel SEQUOIA focal- 
plane array receiver and the DCC autocorrelator, which allowed 
us to observe two different spectral lines simultaneously, either 
in the 85-100 GHz or the 100-1 15 GHz bands. As the two main 
lines of the project, N 2 H + (l-0) and C 18 O(l-0), could not be 
observed simultaneously, the L1495/B213 region was mapped 
twice. In one pass, N2H + (J=l-0) and SO(Jn=32-2i) were ob- 
served simultaneously, and in the other pass, C 18 O(J=l-0) and 
C 17 O(J=l-0) were observed. To achieve high velocity resolu- 
tion, the DCC autocorrelator was configured to provide 1024 
spectral channels with a spacing of 25 kHz, or approximately 
0.07 km s _1 . Rest frequencies are assumed to have the same val- 
ues as in Hacar & Tafalla]( |2011| l. 

All observations were done in on-the-fly mode, covering the 
cloud with a mosaic of 10' x 10' maps referred to a common cen- 
ter at a(J2000) = 4 /, 17 m 47 s .l, 5(72000) = +27°37'18". Position 
switching was used to subtract the sky and receiver contribu- 
tions, and the reference position was located (+2300", -4500") 
with respect to the map center. This position has a C 18 O(J=l-0) 
peak intensity T m b < 0.03 K averaged over the SEQUOIA foot- 
print, as estimated from frequency switched observations. This 
intensity is negligible compared to the typical on-source values, 
which are larger than 1 K. 

During the observations, calibration was done every 10 
minutes, and pointing corrections were determined every 3 
hours from SiO maser observations, finding typical errors 
within 5" rms. According to FCRAO-provided informa- 
tion (http : //www. astro .umass . edu/~fcrao/observer/ 
statusl4m.html), the beam efficiency of the telescope is 
approximately 0.5, and the FWHM beam size depends linearly 
on frequency, having a value of 56" at the lowest (N2H + ) 
frequency and of 47" at the highest (C 17 0) frequency. 

Data reduction consisted in the creation of Nyquist-sampled 
maps with the otftool program (http : / /www. astro .umass . 
edu/~fcrao/library/manuals/otfmanual .html), and re- 
sulted in data sets of more than 35,000 spectra for each of 
our four program lines. These spectra were later converted into 
the GILDAS/CLASS format (http : //www . iram . f r/IRA MFR/ 
GILDAS) for second-degree baseline subtraction and a spatial 
convolution with a gaussian to eliminate residual noise. The final 
resolution of all the FCRAO data, taking into account the off-line 
convolution, is 60", and the grid spacing is 30". With this sam- 
pling, the size of each molecular dataset is approximately 23,000 
spectra. 

Additional line observations of L1495/B213 were made us- 
ing the IRAM 30m telescope. These complementary data con- 
sisted of simultaneous N2H + (l-0) and C ls O(2-l) spectra ob- 
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served toward a selected group of positions to clarify the ve- 
locity structure of the gas. The observations were done in 2010 
April and 2011 December using the EMIR heterodyne receiver 
in frequency switching mode. The backend was the VESPA 
autocorrelator configured to provide a velocity resolution of 
0.06 km s _1 . Sky calibration was carried out every 15 min- 
utes, and pointing corrections every 2 hours. Conversion to the 
mean brightness temperature scale was done using the facility- 
provided telescope efficiencies. The angular resolution of the ob- 
servations was approximately 25" for N2H + (l-0) and 12" for 
C 18 0(2-l). 

Dust continuum observations of selected parts of 
L1495/B213 were carried out at 1200 /mi with the MAMBO-2 
array on the IRAM 30m telescope during several runs between 
2003 December and 2009 November. The observations were 
done in on-the-fly mode with a scanning speed of 8" s , a 
wobbler period of 0.5 s, and wobbler throws between 50" 
and 70". Atmospheric calibration was carried out using data 
from sky dips every 1 .5 hours, and the absolute calibration was 
achieved using observations of CRL618. The data were reduced 
with the MOPSIC software using a method to recover extended 
emission consisting of 10 iterations with a model source and 
subtraction of correlated noise. 

Additional dust continuum observations of L1495/B213 
were carried out at 870 //m with the LABOCA array on the 
APEX telescope in 2007 November and 2009 July. The observa- 
tion consisted of two on-the-fly maps (2007 observations) and a 
mosaic of seven raster-spiral (quarter) maps ( Siringo et al. 2009 ) 
(2009 observations). Calibration was carried out using sky dips 
and observations of Mars and CRL618 about every 2 hours, and 
the reduction used the BOA software together with the method to 
recover extended emission presented by Belloche et al. ( 201 l| l. 
An artistic rendition of these observations can be found in the 
ESO photo release esol209 (http://www.eso.org/public/ 
news/esol209/). 



3. Large-scale cloud properties from the integrated 
maps: evidence for sequential fragmentation 

Figs. [T] and |2]present integrated-intensity maps of L1495/B213 
in C'^O, SO, and N2H/. The different species provide comple- 
mentary information on the gas conditions due to a combination 
of excitation and abundance effects. C 18 is the most sensitive 
tracer of the low-density material in the cloud thanks to its lower 
critical density and relatively large abundance, although it dis- 
appears from the gas phase in dense, chemically evolved regions 
due to freeze out onto dust g rains (Cas elli et al.|1999[|Bergin et] 
al.|2002 Tafalla et al.|2002) >. SO, on the other hand, presents a 



south east, and contains several elongated regions of enhanced 
intensity. Finally, it forks into two branches of relatively diffuse 
emission towards the east and the south-east. 

Overall, the distribution of C 18 emission matches the distri- 



higher critical density that makes it more selective of dense gas 
Its abundance, however, is especially sensitive to the evolution- 
ary state of the material, and is higher during the earliest phases 
of gas contraction (Tafa lla et al.||2006t |Hacar & Tafalla||201 1 1. 
Finally, N2H + is a tracer of dense, chemically evolved gas due 
to its abundance enhancement when CO freezes out, and thus 



highlights the population of dense cores in the cloud (see |Bergin 
& Tafalla 2007 for a review on core chemistry). 

As can be seen in Fig.[T] the C 18 emission (top panel) traces 
a curved filamentary cloud over approximately 4°, or about 10 pc 
for our assumed distance of 140 pc ( Elias|19 78). This emission 
has a bright condensation towards the northern end that coin- 
cides with the location of the LI 495 dark cloud and with one of 



bution of extinction derived by Schmalzl et al. (201011 from NIR 



data with an angular resolution similar to that of our FCRAO 
observations. It also matches approximately the LABOCA and 
MAMBO dust continuum maps presented in the following sec- 
tion. There are however significant departures between the C 18 
emission and both the dust extinction and emission maps. These 
departures coincide with the location of the dense cores that 
dominate the N2FL emission, and they can be seen by comparing 
the maps of C 18 and N 2 H + . Between Aa = 1000" and 3000", 
for example, there is a chain of several cores that is bright both 
in N2I-L and the dust continuum (next section), but that presents 
only weak and diffuse C ls O emission in the map of Fig.[T] This 
anticorrelation between the C I8 and N2FL emissions is typical 
of regions of core formation, and indicates that CO has depleted 
from the de nse gas ( |Caselli et al.[T999"l|Bergin et al.|2002||TafaiIa| 
et al.|2002| . While CO depletion is common in the inner part of 



the main groups of pre-main sequence stars in Taurus (Kenyon et 
|al.|2008| l. From L1495, the emission extends first south and then 



dense cores, the region with the core chain in B213 is unusual in 
presenting depressed CO emission over a scale more than 0.5 pc 
in length. 

While the large-scale geometry of L1495/B213 suggests that 
the cloud is a single entity with a common physical origin, the 
internal structure of the gas presents strong evidence for frag- 
mentation. This can be seen in the maps of Fig. [T] and it was 
already noticed by Barnard from the study of the optical images, 
which lead to his sub-division of the cloud into distinct conden- 
sations (Barnard 1927). To better refer to the different parts of 
the cloud in our study, we have divided the mapped area into re- 
gions, and we have labeled these regions following the original 
notation by Barnard as closely as possible. Barnard's descrip- 
tion, however, is only qualitative, and some parts of the cloud are 
not assigned to any of his condensations, so the correspondence 
between our regions and Barnard's regions should be considered 
only approximate (see also Schm alzl et al.|20lT)| for a similar at- 
tempt to sub-divide L1495/B213). The resulting cloud division 
is shown in the bottom panel of Fig. [2] using color-coded boxes 
superposed to the lowest C 18 O(l-0) contour from in Fig.[T] 

Barnard's division of the cloud was purely based on the ap- 
pearance of the optical images, but is correlated with a number 
of independent cloud properties. One of these properties is the 
stellar population. In Fig. |2j we use star symbols to represent 
the stellar objects classified as "previously identified" or "most 
likely" (A+ rank) Taurus members by the Spitzer telescope sur- 
vey of |Rebull et al. ( 2010) >. Identifying the full stellar population 
of this Taurus region is still subject to debate, especially in the 
low luminosity end of the distribution (e.g., |Palau et aL|2012| l, 
and for this reason we have required conservatively that our can- 
didates have an A+ rank in the notation of Rebull et aL] ( |20T0) l. 
To distinguish between the different stellar classes, we have used 
in Fig.|2]solid star symbols for the youngest objects (Class I and 
Flat, see Rebul fet al.|2010| l, and open star symbols for the more 
evolved objects (Class II and III). 

As it can be seen, the stellar population of the cloud is not 
distributed uniformly over the different Barnard regions. Most 
YSOs are located in regions B7, B213, and B217, which together 
contain 51% of the cloud mass but have 88% of the stars (i.e., 
the other 49% of mass only contains 12% of the stars, see the 
Table [TJ. Even among these three regions, there are significant 
differences in the type of associated stellar objects. B7 and B217 
contain twice as many evolved objects as young ones, while the 
B213 region contains the same amount of young and evolved 
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Aa (arcsec) 

Fig. 1. Integrated intensity maps of L1495/B213 in C 18 O(l-0) (top) and 80(^=32-20 (bottom). The offsets are referred to the 
FCRAO map center at a(72000) = 4' , 17'"47 s :i, 5(72000) = +27°37'18", and the maps have been convolved to a resolution of 75" 
to enhance sensitivity. First contour and contour interval are 0.5 K km s _1 for C I8 O(l-0) and 0.4 K km s for SO(32-2i). 
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N 2 H (1-0) 




10 8000 6000 4000 2000 

Aa (arcsec) 

Fig. 2. Top: Integrated intensity map of L1495/B213 in N9H + (l-0) adding all hyperfine components. The white contour corresponds 
to 0.5 K km s _1 . Central position and offsets are as in Fig!JI] First contour and contour interval are 0.5 K km s _1 . Bottom: Schematic 
view of the L1495/B213 complex indicating the regions defined in this paper, which match approximately those of Barn ard] ( |1927| l. 
The black solid line represents the lowest C 18 O(l-0) contour in Fig. n] and delineates the boundary of the cloud. The red lines 
represent the N2H + (l-0) emission, which traces the dense cores, and the red labels identify the cores described in Sect. |4]and 
summarized in Table [2] The star symbols correspond to stellar objects from the survey of |Rebull et al. ( 20 1 Q) > . Solid symbols 
represent the youngest objects (Class I and Flat) and open symbols represent evolved objects (Class II and Class III). 
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objects. These differences suggest that star formation has not oc- 
curred simultaneously in all the regions, but that different parts 
of the cloud have been forming stars at different rates. 

Another property that correlates with the sub-division of 
the cloud in regions is the chemical composition of the gas. 
Comparing the intensities of C ls O, SO, and N2H + in Figs. [I] 
and|2j it is clear that the emission ratio for the different species 
changes between regions, even when taking into account the 
different excitation requirements of the tracers. The most strik- 
ing contrast between regions occurs in the contiguous B21 1 and 
B213 regions. As mentioned before, the B213 region presents 
relatively weak C 18 and SO emission while is associated with 
bright N2H + . This is indicative of strong CO and SO depletion, 
and suggests that the gas in B213 is chemically evolved and has 
condensed into a number of dense cores. These two properties 
are in good agreement with the finding that B213 has the highest 
fraction of young protostars. 

In contrast with B213, the nearby B21 1 region presents very 
bright emission in both C 18 and SO, together with intense 
dust millimeter continuum (see next section). Its N2I-F emission, 
however, is almost undetected by our observations, except for a 
single core (number 7 in the Fig.[2]l that is not associated with the 
main part of the filament (see Fig. [4] below). This combination 
of bright C ls O and SO together with weak N 2 H + indicates that 
the gas in B211 has an unusually young chemical composition. 
Such chemical youth is again in good agreement with the lack 
of stellar population, which has already been used to argue that 



the B211 region is at a very early state of evolution (Goldsmith 
letaLT2008 1 Schmalzl et al.12010). 



While B211 presents the strongest evidence for chemical 
youth, BIO also seems to be relatively un-evolved. This can be 
seen in Figs. [T]and|2j where BIO presents relatively bright SO 
emission, especially when compared to the C ls O-brighter B7 re- 
gion to its north. In contrast with B211, however, BIO is rel- 
atively weak in C I8 and has a number of N2H + -bright dense 
cores, which suggests that some CO depletion has already taken 
place. BIO must therefore be somewhat more evolved than B2 1 1 . 
Interestingly, while BIO has already formed a number of dense 
cores, none of them seems associated to any Class I or Flat ob- 
jects in the catalog of Rebull et al.| ( 2010 >. 

In summary, the evolutionary state of the different Barnard 
regions is as follows. B211 appears to be the least evolved re- 
gion due to its chemical composition, lack of protostars, and 
presence of a single (starless) core. BIO seems to follow in terms 
of evolution, since it has a number of dense cores but only one 
is possibly associated with a YSO B216 appears also to be at a 
similarly young stage, due to its lack of protostars and cores, but 
its more diffuse nature makes its status less clear. Next in evo- 
lution comes B213, that has started recently forming stars and 
is doing so very actively. Follows B218, which has a number of 



starless cores, including the very young L1521E ( |Hirota et al. 



2002 1 Tafalla & Santiago 2004), and has a very young stellar ob- 



ject. Finally, B7 and B217 have a number of Class I/Flat objects, 
but whose stellar population is dominated by more evolved ob- 
jects. Among these two regions, B7 has the largest number of 
Class III objects (7 out of 25 YSOs), so it is likely to be more 
evolved than the rest. 

While the differences in evolutionary stage between the parts 
of the L1495/B213 complex are clear, the assignment of relative 
ages is more uncertain. Overall, we can estimate that some parts 
of the cloud have been forming stars for at least 1-2 Myr. B7, 
for example, has 7 Class II and Class III objects with ages in 




Aa (c 



3) 



Fig. 3. Expanded view of the B7-B10 region. C 18 O(l-0) emis- 
sion in color and N2H + (l-0) emission in white contours (sum of 
all hyperfine components). First contour and contour interval are 
0.5 K km s _1 . Star symbols as in Fig. [2] and offsets referred to 
the FCRAO map center (see Sect. |2]i. 



age of 2.5 Myr, and B217 contains FV Tau/c, which seems older 
than 1 Myr ( |Kraus & Hillenbrand|2009| l. At the other end of the 
scale, the B211 region, despite being more massive than B213 
or B217, seems to have not yet formed a single star. This means 
that the large-scale cloud, despite its striking appearance as a 
single 10 pc filament, has not fragmented as a single entity, but 
that different parsec-sized regions have followed different star- 
formation histories over the last 1-2 Myr. 

Since the differences in stellar population are correlated with 
differences in the chemical composition of the gas (Table[T]i, and 
since this composition is very sensitive to the amount of time the 
gas has remained dense (e.g., Bergin & Tafalla|2007| l, it seems 
unlikely that regions with little or no star formation, like B21 1, 
have been "waiting" with their current physical conditions for 
more than 1-2 Myr, while regions like B217 were forming stars. 
To appear chemically young, regions like B211 must have con- 
densed more recently from less-dense cloud material than re- 
gions like B217. This suggests that the complex fragmentation 
of the cloud reflects a similarly complex pattern of assemblage of 
the large-scale filament, by which "evolved" regions like B7 and 
B217 condensed and started to form stars first, while other re- 
gions like B21 1 seem to have become dense only very recently. 
As we will see below, this pattern of differential evolution has 
left an imprint in the large-scale pattern of gas velocities. 



the range 1 -6 Myr, according to the estimates by |Bertout et al. 
( |2007| l. B213 is associated with FS Tau B, with an estimated 



4. Dense cores in L1495/B213 

The maps in Figs. [T] and [2] are too large to reveal the details 
of the compact NiH 1 ^ emission and to allow a comparison with 
C 18 0. In Figs. [3|6] we zoom in towards the regions of bright 
N2H + emission and superpose this emission (in white contours) 
with the emissions of C 18 and the dust continuum if available 
(color). 
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Table 1. Properties of the sub-regions in L1495/B213. 



Region 



B7 

BIO 

B211 

B213 

B216 

B217 

B218 



Mass ' 
(Mo) 



Cores <2) YSOs <3) I(SO)/I(C I8 0) (4) Cores/YSOs I+F/II+III <5) 



205 
70 
138 
114 
83 
43 
52 



25 
1 


14 

2 
5 
3 



0.18 
0.42 
0.40 
0.22 
0.28 
0.13 
0.23 



0.12 
5.0 
inf 
0.5 
0.0 
0.0 
1.0 



0.4 
0.0 

1.0 

0.0 
0.3 
0.5 



Rebull et al. 



(1) From C ls O assuming an abundance of 1.7 x 10 7 and T K = 10 K; (2) as determined in Sect. |4j (3) as determined by 
(4) ratio of mean SO(32-2i) and C 18 O(l-0) intensities; (5) ratio of Class I+Flat sources over Class II+Class III sources according to 

|aI~ll |20T0) . 



(2010); 



Rebull et 



As discussed before and clearly seen now in Figs. [3]6J the 
C ls O emission avoids systematically the regions of brightest 
N2H + . This was attributed to CO depletion, a conclusion that 
we confirm from the good match between the N2FF emission 
peaks and the peaks of the dust continuum. This match guaran- 
tees that the N2H + peaks arise from true enhancements of the gas 
column density, and rules out a possible explanation in terms of 
anomalous excitation of N2I-F. Specially striking is the chain of 
four cores shown both in the B21 1 and B213 panels (towards the 
SW in the first one and towards the NE in the second one). This 
0.5 pc-long chain is very prominent both in N2H + in the dust 
continuum, but is practically invisible in C ls O. 



4. 1 . Core selection and physical properties 

As a first step in our analysis of the core population in the 
L1495/B213 complex, we make a census of N2H + condensa- 
tions. Focusing on the N2H + emission guarantees an homoge- 
neous core census, since the continuum maps (the alternative in 
a search for cores) do not cover the full extent of the cloud. It 
however results in us missing a small fraction of N2H + -poor con- 
densations, like L1521E in B218, that are known to be at the ear- 



Table 2. Dense cores in L1495/B213. 



liest stages of evolution ( Hir ota et al.[2 002 ; Tafalla & Santiago 
[20041 . 

To search for N2H + condensations, we have first identified 
all N2H + peaks with intensity exceeding 1.2 K km s _I and fit- 
ted them with 2D gaussians. For that, we have used the fitting 
algorithm of the MOPSIC program (http://www.iram.es/ 
IRAMES/mainWiki/CookbookMopsic), which determines au- 
tomatically the position, size, and intensity of each emission 
peak. In regions where several cores overlap or are very close 
to each other, the fitting process has been applied sequentially. 
In a first step, the brightest peaks were fitted and subtracted from 
the image. Then, a new search for peaks was carried out in the 
residual image, fitting additional gaussians to peaks that still ex- 
ceeded the intensity threshold. This process was repeated until 
no peaks brighter than the threshold remained. 

Table [2] presents the results of our core search. In total, 19 
N2FF cores were identified in the region. This number likely rep- 
resents a lower limit to the number of cores, because if we use an 
intensity threshold lower than the chosen 1.2 K km s _1 , the core- 
finding algorithm identifies a slightly larger number of cores. 
Visual inspection of these weaker cores, however, shows that 
their determination is uncertain due to the presence of weak and 
extended N2H + emission in several parts of the cloud. To guaran- 
tee that each of our selected cores corresponds to a well-defined 



ID 


cy(J2000) 


(5(J2000) 


I(N,H + ) 


D (D 


YSO 


ON' 2 ' 




in s\ 


C ' ") 


(Kkms" 1 ) 


(") 






1 


04 17 43 


28 08 03 


1.6 


140 


N 


5 


2 


04 17 50 


27 56 07 


1.6 


81 


N 




3 


04 17 56 


28 12 23 


1.8 


149 


N 




4 


04 18 04 


28 08 14 


1.4 


114 


N 


8 


5 


04 18 04 


28 22 34 


1.2 


87 


N 


7 


6 


04 18 06 


28 05 41 


2.1 


157 


N 


8 


7 


04 18 10 


27 35 29 


1.8 


135 


N 


9 


8 


04 18 34 


28 27 37 


2.8 


177 


Y< 3 > 


11 


9 


04 18 41 


28 23 22 


1.2 


178 


N 




10 


04 19 37 


27 15 48 


2.1 


139 


N 


13a 


11 


04 19 44 


27 13 36 


2.8 


123 


Y< 4 > 


13b 


12 


04 19 52 


27 1 1 42 


1.9 


100 


N 




13 


04 19 59 


27 10 30 


2.0 


116 


Y< 5 > 


14 


14 


04 20 15 


27 05 59 


1.1 


171 


N 




15 


04 20 59 


27 02 29 


1.8 


258 


Y< 6 > 


16b 


16 


04 21 21 


27 00 09 


2.6 


156 


N 




17 


04 27 54 


26 17 50 


2.6 


137 


N 


26b 


18 


04 28 02 


26 19 32 


1.2 


99 


Y< 7 > 


26b 


19 


04 28 14 


26 20 34 


2.1 


134 


N 





(1) FWHM of the N 2 H + emission (uncorrected for 60" beam); 

(2) Core number in the catalog of |Onishi et al. ( |2002| >; 

(3) 104152+2820; (4) 104166+2706; (5) 104169+2702; (6) 
2MJ042 10795; (7) 104248+2612 



emission peak, we have preferred to use a conservative criterion 
based on a relatively high threshold. Our analysis, therefore, will 
focus on this set of the brightest dense cores of the cloud. 

The results of our core search presented in Table|2]show that 
the L1495/B213 condensations have sizes and N2H intensities 
typical of the Taurus core population (Benson & Myers 1989; 
|Caselli eFa l. 2002). Some of these cores are associated with em- 
bedded Class and Class I objects and therefore must have al- 
ready undergone gravitational collapse, while others are starless 
and likely represent pre-stellar condensations. A number of our 
N?H+ -select ed cores have counterparts in the catalog by Onishi 
|et al.| ( |2002| l, which is based on H 13 CO + observations. Our core 
list, on the other hand, has very little overlap with the larger list 
of extinction-selected cores presented by Schmalzl et al. ( 2010] >. 
As can be seen in the Table 2 of these authors, these extinction- 
selected cores have densities on the order of 10 4 cm 3 , and there- 
fore likely represent a population of condensations at an earlier 
state of contraction than our N2H + -selected cores. 
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< 




2000 1000 

Aa (arcsec) 

Fig. 4. Expanded view of the B211 region and the NW part of 
B213 (to better connect with the following figure). Top: C ls O(l- 
0) and N2H + (l-0) maps as in Fig. [3] Bottom: 850 pm contin- 
uum map from LABOCA observations and N2H + (l-0) emission 
(contours). 



4.2. Core clustering 

As we saw in Fig. [2) the N 2 H + cores in L1495/B213 are not uni- 
formly distributed over the cloud, but seem to form small groups 
of a few objects each. A similar trend for clustering, but in the 
lower density gas, was found by Schmalzl et al.| ( |2010| l using 
a two-point correlation analysis of the distribution of extinction 
peaks. In this section, we study the clustering of the N2H + cores 
using the mean surface density of companions S(0), which is 
equivalent to the two-point correlation function but more often 
used in studies of stellar clustering (e.g., Larson 1995 see also 
Simo n|1997|f or the simple relation between the two). Following 
Simon| | |1997) , we have measured for each core the angular sep- 
aration to all the other cores in the cloud, and we have binned 
the set of angular separations in logarithmic intervals with a step 




Aa (arcsec) 

Fig. 5. Expanded view of the B213 region. Top: C I8 O(l-0) and 
N2H + (l-0) maps as in Fig. [5] Bottom: 1200 pm continuum 
map from MAMBO observations and N2H + (l-0) emission (con- 
tours). 



equal to a factor of 2. To convert this number of separations into 
£(#), we have divided the result by the area of the bin (which 
we have assumed rectangular with a width of 100"), and finally 
we have normalized the value by the total number of cores in the 
cloud. 

The resulting density of core companions in L1495/B213 
is presented in Fig. [7] with solid red squares. To compare it 
with the expectation from a truly uniform distribution of cores, 
we have carried out a set of 100 Monte Carlo simulations. In 
each simulation, 19 different cores have been created with co- 
ordinates randomly distributed over a rectangle with dimensions 
13100" x 100", which is approximately the size of the C 18 0- 
emitting region in Fig. [T] Each of these sets of 19 cores has 
been treated like the original N2H + data set, and has been used 
to calculate a separate mean surface density of core compan- 
ions with the same code used for the data. The mean of the 
100 Monte Carlo experiments is represented in Fig.|7]with open 
black squares together with error bars indicating the rms value 
of the dispersion. 

As can be seen from the figure, the N2H + data and the Monte 
Carlo simulation agree within approximately l-cr in the 4 largest 
bins, indicating that the distribution of cores with separations 
larger than about 700" (0.5 pc) is consistent with being spatially 
uniform. At smaller distances, however, the observed distribu- 
tion of cores deviates from the prediction of the uniform model. 
This is specially striking for second smallest bin, which is cen- 
tered at core separations of 350" (0.25 pc). In this bin, the ob- 
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0.5 pc 




8000 

Aa (arcsec) 



Fig. 6. Expanded view of the B218 region. Top: C 18 O(l-0) and 
N2H + (l-0) maps as in Fig. [3] Bottom: 850 pm continuum map 
from LABOCA observations and N2H + (l-0) emission (con- 
tours). 
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Fig. 7. Mean density of companions for the N2H + cores pre- 
sented in Table [2] The red solid squares represent observations 
and the open black squares with error bars represent the predic- 
tion from a series of 100 Monte Carlo runs (see text). 



served density of core companions exceeds the model predic- 
tion by a factor of three, indicating that our sample has three 
times more cores with separations of 2-3 core diameters than ex- 
pected for a uniform distribution of objects. Whether a similar 
or smaller excess of cores occurs at smaller separations is un- 
clear. The data in the first bin of Fig. [7] presents an excess that 
is not statistically significant, but our experiments with the core- 
finding algorithm show that this bin is prone to incompleteness 
due to our limited angular resolution and to confusion with ex- 
tended emission. Higher quality data are required to analyze the 
behavior of the core distribution at distances of 1-2 core diame- 
ters. 

The excess of cores with separations of about 0.25 pc is 
a strong indication that the conditions of core formation in 
L1495/B213 are not randomly distributed over the cloud. Some 
regions seem to be specially favorable to form cores, and there- 
fore give rise to multiple condensations in close proximity. 
These regions, however, cannot extend for much more than about 
0.5 pc, as the mean surface density of core companions approxi- 
mates a uniform distribution at scales of that size or larger. As we 



will see below, this behavior of the distribution of cores can be 
understood if core formation occurs by the fragmentation of fila- 
mentary structures that have 0.5 pc or so in length, provided that 
the cloud is composed of a large number of filaments, and that 
only some of them have conditions leading to core formation. 
The signature of core formation by fragmentation of 0.5 pc-long 
structures seems therefore imprinted in the distribution of core 
positions. 

5. Velocity structure: evidence for multiple 
components 

Even taking into account the large-scale fragmentation of the 
L1495/B213 complex discussed before, the maps of integrated 
emission and dust continuum give the impression that the mate- 
rial in L1495/B213 consists of a single filamentary structure. The 
molecular spectra, however, reveal a much more complex pic- 
ture, and show that the cloud contains multiple velocity compo- 
nents that often overlap in projection. This can be noticed from 
the inspection of the FCRAO C 18 O(l-0) data, but it is better 
appreciated using a series of high signal-to-noise spectra taken 
using the IRAM 30m telescope with the specific goal of charac- 
terizing the velocity structure of the cloud. Some of these spec- 
tra are shown in Fig. [8] which illustrates the kinematics of the 
gas in the vicinity of of the B213 region, the one with the most 
complex velocity pattern. Previous observations of this region 
had noticed the presenc e of two velocity components separated 
by more than 1 km s' 1 (Heiles & Katz|l976||Crark et al.|l977[ 
IDuvert et al.1[T986l [Oriishi et al.||1996f|Li & Goldsrnith||2012| i. 
Our high sensitivity IRAM 30m data reveal now that the kine- 
matics of the region is more complex, and that the number of 
C ls O velocity peaks is higher and changes rapidly with position. 
Fig.[8]shows how each of two previously known velocity compo- 
nents (near 5.3 and 6.7 km s~', see panel 2), splits in some places 
into additional well-separated components, like those illustrated 
in panel 4. A more extreme splitting occurs in the the vicinity 
of (-100", 200"), where the C 18 spectrum presents 5 separate 
peaks (panel 3). Surprisingly, this region with multiple velocity 
peaks corresponds to a single emission maximum in both the 
continuum and the FCRAO C 18 integrated intensity map. 

The multiple velocity peaks in the C ls O spectra of Fig. [8] 
seem to originate from true velocity components in the gas re- 
sponsible for the emission. Self absorption could also create 
multiple velocity peaks in a spectrum by depressing the inten- 
sity at intermediate velocities, but this effect seems to not occur 
in the C 18 data. The isolated component of N 2 H + (l-0) (F\F = 
01-12), which according to hyperfine-structure fits is optically 
thin in 86% of our positions, always appears if detected at the 
velocity of one of the C ls O components, and not between the 
two, as expected in a case of self absorption (e.g., |Leung|1978 1. 
This is illustrated in panels 1, 2, and 4 of Fig. [8] where the iso- 
lated component of N2H + (l-0) is superposed in red to the multi- 
component C I8 spectrum from the same position. As can be 
seen, the N2H + line is associated with one of the C ls O peaks in a 
way that is inconsistent with a self-absorption origin of the C 18 
velocity peaks (the lack of N 2 H + in the other velocity component 
suggests it has a lower density or a lower N2H + abundance). At 
the extreme position with 5 components in the C 18 spectrum 
(panel 3), no clear emission is seen in the isolated component 
of N2H + , but a nearby position less than 60" away presents an 
N2H + spectrum with two weak components that coincide with 
the two reddest peaks of the two groups of C 18 components. 
This again argues for self absorption not being the cause of the 
multiple peaks seen in the C 18 spectrum. 
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Fig. 8. Central panel: FCRAO integrated intensity maps of the B211-B213 region in C I8 O(l-0) (grey scale and black contours) 
and N2H + (l-0) (red contours). Left and right panels: Spectra from selected positions illustrating the complex velocity structure of 
the emission. The blue spectra correspond to C I8 0(2-1) and the red spectra correspond to the isolated {F\F = 01-12) component of 
N2H + (l-0) (multiplied by 3) observed simultaneously with the IRAM 30m telescope. Note the presence of 5 different C 18 peaks 
in panel number 3. 



The multiplicity of velocity components in the C ls O spectra 
presents a significant challenge for representing and analyzing 
the emission. Velocity structure in a cloud is commonly charac- 
terized using channel maps, in which the emission is integrated 
in narrow ranges of velocity that can be selected to highlight par- 
ticular gas components. For the L1495/B213 complex, however, 
this approach fails due to the large number of components, their 
partial overlap in velocity (as shown by the spectra in Fig. [8j, 
and the fact that each component changes slightly its velocity 
with position. A better approach to study the velocity field is 
to search for the different velocity components directly in the 
spectra. Fig.|8]shows that the velocity components are relatively 
symmetric when found in isolation or detected in N2H + , and this 
suggests that fitting multiple gaussians to the emission can deter- 
mine the main properties of each component, notably its average 
velocity and its velocity dispersion. Fitting manually an indeter- 
minate number of gaussians to each spectrum in a dataset that 
contains tens of thousands of spectra is highly impractical, so we 
have developed a semi-automatic procedure. A number of tests 
suggested that the following steps provide a reasonable balance 
between automation and supervision. 



Divide the mapped area in fields of 150" x 150", each con- 
taining 25 individual spectra. 

For each field, discard all spectra with no channel brighter 
than 3 times the typical rms (w 0.3 K). If less than 3 spec- 
tra are left in a given field, the field is discarded from the 
analysis. 



Average all accepted spectra in a field to obtain a high S/N 
template, which is inspected visually to determine the num- 
ber of gaussian components needed to fit it (in cases of doubt, 
the individual spectra were inspected). To be considered real, 
a component needs to lie at least 0.25 km s _] (3 channels) 
apart from the other components in the spectrum. 
Fit each spectrum in a field using the number of gaussian 
components determined in the previous step. This is done us- 
ing the MINIMIZE command of CLASS. The fit starts with an 
initial guess based on the average spectrum of the field, but 
the properties of each gaussian component (intensity, cen- 
tral velocity, and dispersion) are left unconstrained and fitted 
automatically. 

Discard those fits whose intensity is lower than 3 times the 
rms in the spectrum. 



With this procedure, it was possible to process reasonably 
quickly the approximately 23,000 spectra of the FCRAO 
C 18 O(l-0) map and fit a total of more than 1 1,000 gaussian com- 
ponents with SNR equal to or larger than 3. For the N2H + (1- 
0) data, the analysis was simpler, since no FCRAO spectrum 
showed two velocity components with high enough signal to 
noise to deserve a multi-gaussian fit. We therefore applied a 
single-component hyperfine-structure fit (using the HFS option 
in CLASS) and again selected as real those fits exceeding 3 times 
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Fig. 9. Position-position-velocity cube showing the line center 
velocities derived from the gaussian fit to the C 18 O(l-0) and 
N2H + (l-0) spectra towards the B211-B213 region presented 
in Fig. [8] Note the presence of correlated structures. Semi- 
transparent points have intensities between 3 and 6 times rms 
and opaque points have intensities larger than or equal to 6 rms. 



the rms in the spectrum. Almost 400 N2H + spectra were fitted 
this wa>Q 

As each gaussian component is characterized by a position in 
the sky and a center velocity (plus a width), the most convenient 
way to explore the results of the multi -component gaussian fit to 
the data is by plotting the component line center velocities in a 
cube of position-position-velocity (PPV) space. Such PPV space 
represents a subset of the gas phase space and allows visualizing 
structures spatially confused in the plane of the sky but separated 
in velocity. Fig. [9] shows one such a PPV cube that covers ap- 
proximately the same region presented in the integrated map of 
Fig. [8] and contains the result of our gaussian fitting. Although 
undoubtedly complex, the distribution of points in PPV space 
shows clear signs of arising from an organized structure. Most 
points cluster in elongated groups that lie at different "heights" 
(velocities) and present smooth and often oscillatory patterns. 
These groups are associated with the different velocity compo- 
nents seen in the spectra, and their clean separation in the cube 
shows how working in PPV space can help disentangle the com- 
plex gas kinematics of the cloud. 

The next step in our analysis requires a procedure to identify 
and isolate the different components in the PPV cubes. Visually 
exploring those cubes, it is possible to identify a number of indi- 
vidual groups of points that form coherent structures and likely 
correspond to physically distinct cloud components. It seems 
however preferable to have an objective algorithm that not only 
automatizes the selection, but also provides a quantitative mea- 
surement of the degree of correlation between points in PPV 
space. 

6. Friends In Velocity (FIVe): an algorithm to 
identify velocity components 

The problem of identifying velocity structures in a PPV cube 
like that of Fig. [9] has a number of similarities with the prob- 

1 The fits results for all the components with SNR > 3 identified in 
our C 18 (1-0) and N 2 H + (1-0) FCRAO spectra are available online at 
the CDS. 



lem encountered in cosmology of identifying groups of galaxies 
in a redshift survey. Redshift surveys provide catalogs in which 
each galaxy is characterized by two coordinates and a redshift 
(equivalent to our PPV data), and a search for galaxy groups 
involves identifying those sub-sets of galaxies that are strongly 
correlated in both position and redshift. One of the first and sim- 
plest algorithms to identify galaxy groups in redshift surveys is 
the so-called friends-of-friends (FoF) method, initially presented 
by Huchra & Geller| ( |1982[ ) and still widely used (e.g., |Berlind| 



et al. 2006). This FoF algorithm starts by selecting a random 
galaxy in the catalog and searching the rest of the catalog for 
related galaxies ("friends") using pre-determined thresholds in 
both position and redshift. Each of the galaxies identified as a 
friend is added to the galaxy group, and a new search for friends 
is made this time centered on the friend galaxies. The new finds 
(friends of friends) are added to the group, and the process is it- 
erated until no more friends are found, at which point the group 
is considered complete. When this happens, a new unassigned 
galaxy is selected from the catalog and a new search for friends 
is started to identify a new galaxy group. The final product of 
the algorithm is a catalog of galaxy groups and a catalog of iso- 



lated galaxies not belonging to any group (see Fig. 1 in |Huchra 
|& Geller|198"2"l for a flow chart). 

Applying the FoF technique to the line center velocity of the 
C 18 components follows the spirit of the cosmology implemen- 
tation, but requires a number of adjustments due to the different 
nature of the emitting material. Each galaxy in a redshift survey 
is a discrete entity, physically separated from the rest (excluding 
mergers), and therefore well represented by a single point in PPV 
space. A velocity component in a C ls O spectrum, on the other 
hand, represents a parcel of a fluid that extends over a large re- 
gion of space and has no sharp boundaries. The C 18 emission, 
in addition, has been mapped with Nyquist sampling and the data 
points form a dense grid, in contrast with the sparse sampling 
of a redshift survey data. These peculiarities of the C I8 emis- 
sion require modifying the standard FoF approach, in particular 
to avoid a weakness of the method when applied to fluids and 
which we will refer to as "fragility". This fragility of the FoF 
method means that a single bad point (e.g., from a poor gaus- 
sian fit) can potentially create an artificial bridge between two 
components that are otherwise disconnected in PPV space, and 
that represent physically independent entities. While this is not 
likely to occur in redshift data because of their sparse nature, it 
is a real problem when dealing with the Nyquist-sampled and 
diffuse C 18 emission. This fragility can lead, for example, to 
an absurd situation in which two distinct velocity components 
in a spectrum become connected by the FoF algorithm because 
each one is linked through a number of friends to a distant and 
weak point with an intermediate velocity, either due to a poor fit 
or because of a more complex velocity pattern. 

To mitigate the fragility of the FoF algorithm, we need to 
combine the proximity concept of the standard FoF implementa- 
tion with an additional constraint. One possibility is to use the in- 
tensity of the emission as an additional criterion to identify com- 
ponents. This approach seems the most natural extension of FoF 
when dealing with diffuse emission, since intensity itself (which 
is proportional to column density) is the criterion already used 
to define cloud components in velocity-insensitive datasets like 
continuum or integrated-intensity maps. By adding an intensity 
consideration to the FoF algorithm, our goal is to make it capa- 
ble of recognizing that structures like two long ridges separated 
by a deep emission valley represent different cloud components, 
even if the velocity of the gas in both ridges differs by less than 
the FoF-imposed threshold. Adding intensity weighting to the 
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a) Step 1: High-density points 
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b) Step 2: Friends-of-Friends algorithm and identification of 
velocity-coherent structures 
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c) Step 3: Association 
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Fig. 10. Schematic view of the 3 steps in the modified FoF 
method used to identify components in PPV space. To simplify 
the view, two projections of the PPV cube are presented. The left 
panels show the PP plane (integrated map) and the right panels 
show a projection on a PV plane. The black points represent all 
data with SNR > 3. The green points in the first row of panels 
are those selected as significant to use in the first application of 
the FoF algorithm. The red and blue colors in the lower panels 
indicate positions belonging to each of the two structures iden- 
tified by FoF as connected. The colored lines indicate the indi- 
vidual connection between points. Note how points are gradualy 
assigned to the two velocity structures. The boxes illustrate the 
regions used to calculate the number of neighbor points and the 
velocity gradients for the FoF algorithm. 



FoF algorithm is however a complex task. A number of geomet- 
rical and intensity considerations are needed, and the sophisti- 
cation of models used to analyze intensity-only continuum data 
from the Herschel Space Observatory testifies the needed level 
of complexity (e.g., |Men'shchikov et al.|2012) . To analyze our 
data, we have chosen a simplified method that merges the con- 
cepts of FoF and intensity weighting, and that has been shown 
to work well in the PPV data set produced from our FCRAO ob- 
servations. The main idea behind the method is that the emission 
intensity can be used to apply the FoF algorithm selectively in a 
top-down manner, starting from the brightest points in the map 
and working the way downwards in sensitivity. A cutoff that se- 
lects the brightest 50% points in the dataset is used to define the 
independent components in PPV space. This approach, which is 
illustrated in Fig. 10 with an idealized case, consists of the fol- 
lowing three steps. 
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b) Step 2: FoF algorithm and identification of Velocity-Coherent Structures 
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Fig. 11. Analysis of the PPV cube of Fig^\ using the FIVe al- 
gorithm. Top: Seed points selected using the threshold values 
described in the text (color coded green). Middle: Points iden- 
tified after a friends -of-friends search that starts with the seed 
points from the previous step. The color coding illustrates the 
different cloud components. Bottom: Final assignment of points 
to the different cloud components (colored points). The black 
points represent unassigned points. 



We first identify the most significant points in the PPV cube 
that can serve as seeds for the FoF algorithm. We do so by 
finding positions that either have a SNR > 3 in N2H + , or 
that have both SNR > 6 in C 18 C0 and more than half of 
the 8 nearest neighbors in the 30" grid with SNR > 6 and 



2 The intensity of the C 18 O(l-0) emission in positions with bright 
N2H + (l-0) lines (SNR > 3) has been multiplied by 2 to compensate for 
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a difference in velocity equivalent to a gradient lower than 
3 km s _1 pc -1 . (A 3 km s _1 pc -1 threshold guarantees that 
no two points in the 30" differ in velocity by more than the 
typical non-thermal dispersion of the C 18 O(l-0) line, which 
is on the order of the sound speed, or 0.19 km s .) These 
criteria were found to select the w 35% of points in the sam- 
ple that lie in regions of homogeneous properties from which 
the different cloud components can be identified. 

2. Using the points selected in the first step, we run a FoF search 
using the C 18 data with a box of 120" and the previous 
3 km s pc~' velocity gradient threshold. The choice of the 
box size, which covers two layers of the 30" grid around the 
point, allows the algorithm to explore larger distances than 
simply the next neighbors, and provides an insurance against 
small-scale drops in the intensity of the emission that could 
otherwise fragment the component. Once the FoF search has 
been completed, the cloud velocity components are defined 
as those groups of friends that are isolated from the rest and 
that contain a minimum number of 8 members. 

3. The final step in the analysis consists of assigning the points 
with few neighboring friends and/or SNR less than 6 to 
the different cloud components defined in the previous step. 
For this, we relax the previous neighbor and SNR thresh- 
olds, and test each component point for new friends using 
again the 3 km s _1 pc 1 velocity gradient threshold. In line 
with the intensity-weighting scheme discussed before, this 
new search for friends is applied sequentially, advancing 
each cloud component by one step. When no new friends 
are found, the analysis is stopped and the remaining points 
(typically < 20%) are considered unassigned. Most of these 
points are low S/N data or points that belong to what appear 
to be separate cloud regions that did not have enough points 
in step 2 to qualify as a bona-fide component according to 
our definition. 

All the above steps are carried out automatically using a 
dedicated algorithm called FIVe (Friends In Velocity) and pro- 
gramed using the R language (http : //www . r-pro j ect . org/). 
More information on the algorithm and its implementation can 
be found in Hacar (2012). A series of tests to determine its suc- 
cess rate and to illustrate the choice of internal parameters are 
presented in Appendix [A] 

7. The nature of the cloud components 

7.1. Components as filaments 

While Fig. [10] showed an idealized view of our FoF analysis, 
Fig. [TT] presents the true result of applying the FIVe algorithm 
to the the PPV cube of data from the region previously shown in 
Fig. [9] There may be some ambiguity in the assignment of the 
weakest points to the different components, but the algorithm 
seems to identify as separate entities most structures in PPV 
space that visually appear as distinct components in the cube. 
Over the whole cloud, the FIVe algorithm identifies 35 differ- 
ent components, 17 of them in the region shown in Fig. 1 1 The 



main physical properties of these components are summarized 
in Table. [3]and discussed in the rest of this section. 

The most noticeable property of the L1495/B213 compo- 
nents is their filamentary geometry. This can be seen directly in 
the PPV cubes, where many of these components appear as long 
and twisted structures. It is however better appreciated using 

the drop in C 18 intensity due to freeze out, which is typically a factor 
of 2 ([Tafalla et al.|2004|>. 



standard position-position maps. Using these maps, we have de- 
termined for each component a principal axis by connecting with 
straight lines the emission centroids of 90" -long fragments, as 
explained with more detail in Appendix [B] These principal axes 
represent the backbones of the components, and they are shown 



in Fig. 12 superposed to the 250 /jm dust continuum Herschel- 
SPIRE archive from the Gould Belt Survey ( ] Andre et al.|2010 
Palmeiri m et al.|2013) >. 



The most striking aspect of Fig. 12 is the the intricate pattern 



of intertwining filamentary components, especially towards the 
B21 1-B213 regions shown with detail in the figure inset. A mul- 
tiplicity of components was of course expected from the multiple 
velocity peaks seen in the spectra of Fig. [8] The map of princi- 
pal axes, however, shows that the components, despite their very 
different velocities, belong to a rather organized structure. What 
at first sight seemed like a single 10 pc-long filamentary cloud 
becomes now a complex network of braided filaments. 

Another remarkable aspect of Fig. 12 is how well the prin- 
cipal axes of the velocity components follow the 250 fim dust 
emission mapped with SPIRE. This is specially noticeable in the 
B2 11-213 region, where most of the components are located. It 
suggests that the cloud components traced by the combination 
of C 18 and N2H + emission correspond to the same material 
traced by the SPIRE observations of the dust emission. Indeed, 
the mass estimated by Palmeirim et al. ( |2013) > for the B211- 
213 region matches closely the mass estimated by our analysis 
(Table [TJ. The additional kinematic information provided by our 
line data show that the continuum-emitting material consists of 
multiple velocity components, and this multiplicity, as discussed 
in the next section, has important consequences in the estimate 
of the mass per unit length. 

7.2. Statistics of filament properties 

We start the analysis of the filaments using a statistical approach. 
Fig. 13 presents a series of histograms showing how the most rel- 



evant properties are distributed among the 35 filaments identified 
by the FIVe algorithm. In agreement with the previous discus- 
sion, the histogram of aspect ratios (top left) shows a distribution 
significantly shifted from 1, with 75% of the components having 
an aspect ratio larger than 3. This large value of the aspect ratio 
provides a justification for our using the word "filament" to refer 
to most of the observed components. 

The statistics of linear sizes is illustrated in the top-right his- 
togram of Fig. 13 which shows that despite a significant tail of 



large values, 60% of the cloud components lie in a relatively 
narrow peak between 0.2-0.6 pc in sizeP] 

A related quantity is the mass per unit length, which is pre- 
sented in the middle-left histogram (Fig. [13};). This parameter 
has been calculated directly fro m the C 18 Q emission assuming 
the standard C 18 abundance of Frerking et al. ( 19821. To com- 
pensate for CO depletion, the contribution from positions with 
bright N2H + emission (SNR > 3) has been supplemented with a 
mass estimated using the N2H + emission and assuming standard 
ex citation and abundance values (Cas elli et al.| [2002; Taf alla et] 
|aLl|200 4). These N2H + -derived contributions increase the fila- 
ment mass typically by 25% and never reach 50% of the C 18 0- 
derived mass. As the histogram shows, the mean mass per unit 



3 Some of the longest components may result from the artificial 
merging of smaller units. Indeed, the longest component (number 20) 
often splits into two separate components when we run FIVe using 
slightly different thresholds for the FoF algorithm. Higher sensitivity 
data are required to solve this ambiguity. 
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WJ2000 



Fig. 12. Location of the cloud components of L1495/B213 identified with the FIVe algorithm. Each component is represented by 
its central axis and has been color-coded for easier identification. The background grey-scale image is a SPIRE 250pm continuum 
map from the Gould Belt Survey (An dre et al.|20 10) that has been blanked out to match the coverage of our FCRAO data. 



length of the components is about 15 M Q pc 1 . Thi s value is sig- 



nificantly lower than the * 54 M Q pc derived by |Palmeirim et 
[aT7| ( |2Q13] > from their Herschel continuum data, but as mentioned 
before, this difference results from the fact that the continuum 
analysis associates all the mass to a single cloud component, 
while our line analysis shows that the mass is distributed among 
distinct velocity components. 

Interestingly, the derived mean mass per unit length of 
15 M Q pc~' is very close to the equilibrium value for an isother- 



mal cylinder in pressure equilibrium at 10 K ( Stodolkiewicz 
1963 Ostriker|1964| l. This property, and both the length and as 



pect ratio discussed before, make the L1495/B213 components 
very similar to the velocity-coherent filaments of the nearby 
L1517 cloud, also studied using C 18 data. These L1517 fila- 
ments had aspect ratios of approximately 4, typical lengths of 
0.5 pc, and mass per unit length in agreement with the predic- 
tion for a 10 K isothermal cylinder ( |Hacar & Tafalla|201 1) . 

A notable difference with the L1517 filaments is the larger 
velocity dispersion of the gas in the L1495/B213 components. 
The middle-right histogram in Fig. 13 shows the distribution 



of the non-thermal velocity dispersion normalized to the sound 
speed at 10 K, as calculated by taking the unweighted mean of 
the non thermal linewdiths in the individual gaussian fits to the 
spectra (the dispersion of the individual line center velocities 
presents a similar behavior). While the gas in the L1517 fila- 
ments was overwhelmingly subsonic (98% of the points) with a 
typical cr NT /c s = 0.54 ± 0.19, the gas in L1495/B213 presents 
an approximately sonic dispersion of (Tnt/c s = 1.0 ± 0.2, al- 



though there is a ~ 15% minority of subsonic, L1517-like fil- 
aments. This larger value of the C 18 non-thermal motions in 
L1495/B213 does not arise only from the action of the embed- 
ded protostars (whose effect is noticeable but local), and indi- 
cates that the gas in most L1495/B213 filaments is more tur- 
bulent, or has larger internal velocity gradients, than the gas in 
L1517. Even so, the estimated ctntIc s values indicate that the 
internal motions in the L1495/B213 filaments are at most mildly 
transonic, and therefore inconsistent with a possible picture of 
supersonic turbulence. 

Probably related to this larger C 18 dispersion is an also 
larger dispersion of the non-thermal component of N2H + . The 
mean cr NT /c s value for this species is 0.61 +0.17, compared to 
0.36 ± 0.09 in L1517. Apart from this higher velocity disper- 
sion, the behavior of the gas in the L1495/B213 filaments is very 
similar to L1517, as we discuss in the following section. 

The final two histograms present the distribution of velocity 
gradients measured along the long axis of each filament (bot- 
tom panels of Fig.[T3]>. The left histogram shows the distribution 
of global velocity gradients, which are determined by fitting the 
full set of LSR velocity values along the long axis with a single 
linear gradient. The right histogram (Fig. \13\ ) shows the distri- 
bution of local velocity gradients, which are determined using a 
similar method, but this time splitting the velocity data in 0. 1 pc 
fragments and then taking the average of the fits. As can be seen, 
the mean value of the global gradients is close to 0.5 km s -1 pc _1 , 
and the mean value of the local gradients is about 1 km s _1 pc -1 . 
The larger value of the local gradients results from the presence 
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Table 3. Main properties of the L1495/B213 cloud components. 



Fil ID 


Aa" 


AS" 


Mass 


L* 


M,„, £ 


<v,.„>'' 


o-V, s ,-/c/ 


(<T NT )/c/ 


|VV, S ,-| S 


Cores* 




(arcsec) 


(arcsec) 


(M Q ) 


(pc) 


(M G pc-') 


(km s" 1 ) 






(km s _1 pc -1 ) 




01 


-654 


3204 


1.5 


0.1 


11.3 


5.9 


0.3 


0.82 


1.11 





02 


-420 


2545 


30.9 


1.0 


32.4 


7.0 


1.1 


0.99 


0.23 





03 


-181 


350 


1.8 


0.4 


3.9 


5.7 


0.3 


0.86 


0.02 





04 


-136 


228 


1.2 


0.2 


5.4 


5.3 


0.3 


0.81 


0.62 





05 


-130 


316 


4.0 


0.4 


9.1 


4.9 


0.4 


0.73 


0.16 




06 


-99 


1911 


5.3 


0.5 


11.4 


6.0 


1.0 


1.41 


1.29 


1 


07 


-50 


141 


0.7 


0.2 


4.6 


7.0 


0.4 


0.69 


0.91 





08 


-35 


1573 


2.5 


0.3 


9.3 


6.7 


0.5 


1.37 


0.02 





09 


159 


-237 


4.3 


0.5 


8.7 


5.0 


0.7 


0.83 


0.18 





10 


183 


2055 


15.5 


0.8 


20.4 


6.8 


0.7 


1.39 


0.34 


3,4,6 


11 


305 


-44 


6.9 


0.5 


14.7 


6.7 


0.7 


1.27 


0.39 


7 


12 


433 


-502 


29.3 


1.4 


21.8 


5.6 


1.1 


1.25 


0.16 




13 


462 


-981 


1.2 


0.3 


4.3 


6.5 


0.5 


0.84 


1.10 




14 


633 


-650 


5.6 


0.5 


10.3 


5.0 


1.0 


0.97 


0.04 





15 


675 


2712 


95.9 


1.3 


71.4 


7.2 


0.8 


0.97 


0.01 


8.9 


16 


993 


-1156 


1.6 


0.4 


4.5 


4.8 


0.8 


0.81 


1.28 




17 


1114 


3092 


4.3 


0.3 


13.1 


7.2 


0.5 


0.84 


0.38 





18 


1123 


-1255 


15.7 


0.9 


17.0 


5.6 


1.3 


1.01 


0.84 





19 


1483 


-1766 


18.9 


1.1 


17.1 


5.9 


0.8 


1.33 


0.06 





20 


2259 


-2015 


48.3 


2.5 


19.7 


6.6 


1.2 


1.15 


0.11 


10-16 


21 


2637 


-1960 


5.7 


0.6 


9.1 


5.9 


0.7 


1.07 


0.27 





22 


3347 


-2530 


3.1 


0.6 


5.5 


7.0 


0.6 


0.87 


0.34 





23 


3472 


-2565 


4.6 


0.5 


9.4 


6.1 


0.7 


1.05 


0.14 





24 


3630 


-2999 


0.7 


0.3 


2.7 


6.3 


0.4 


0.58 


0.92 





25 


4106 


-2975 


1.0 


0.2 


4.6 


6.4 


0.3 


1.09 


0.35 


— 


26 


4371 


-3255 


1.4 


0.2 


7.7 


7.0 


0.6 


0.77 


1 .55 




27 


4740 


-3248 


1.3 


0.3 


5.0 


6.7 


0.4 


0.76 


0.82 




28 


5227 


-3665 


43.7 


1.8 


24.1 


6.5 


1.0 


1.03 


0.07 




29 


7206 


-4938 


9.8 


0.7 


14.1 


6.8 


0.5 


0.90 


0.34 




30 


7369 


-5465 


3.9 


0.4 


9.1 


6.5 


0.5 


1.17 


0.35 




31 


7369 


-4301 


4.9 


1.0 


5.0 


7.1 


0.7 


0.63 


0.51 




32 


7979 


-4642 


8.7 


0.5 


17.8 


7.1 


0.7 


0.98 


0.12 


17,18 


33 


8257 


-4251 


0.6 


0.2 


3.1 


6.4 


0.4 


0.70 


1.32 




34 


8578 


-4493 


13.3 


1.0 


13.4 


6.9 


1.1 


0.90 


0.65 


19 


35 


9349 


-4960 


5.5 


0.4 


15.8 


6.7 


0.9 


1.11 


1.66 





Notes. ((!) Intensity-weighted centroid measured in offset from the FCRAO map center; length of component determined as explained in 
Appendix|B] (c) mass per unit length; (rf) mean velocity of emission; (e) dispersion of the line center velocity in sound speed units; ( -° mean of the 
nonthermal velocity dispersion in sound speed units; (g) mean gradient of the line center velocity; (/,) associated cores labeled as in Tabled 



of velocity oscillations inside the components, which average 
out when taking a global gradient. A similar behavior was seen 
inL1517. 

7.3. Core formation in the L1495/B213 filaments 

In section [4] we identified 19 dense cores in the L1495/B213 
complex based on their N2H + emission. Now we study the re- 
lation between these dense cores and the 35 filamentary compo- 
nents identified from the C I8 data. Table [3] shows the assign- 
ment of the cores to the different filaments based on the matching 
of the N2H + and C ls O emission in both position and velocity. 
Cores 2 and 5 were not assigned to any filament because their 
surrounding C ls O emission did not meet the minimum number 
of points required for filament definition in the FIVe analysis. 
This lack of association seems a sensitivity problem and not an 
indication that these cores are peculiar, since a FIVe analysis 
using a SNR threshold of 2.5 (instead of the standard 3 value) 
identifies C 18 filaments associated to these cores (but also suf- 
fers from confusion in different regions). 

A notable aspect of the assignment of cores to filaments is 
that only 5 out of 19 cores can be classified as single, in the 
sense that they belong to filaments containing no other core (to 
be conservative, unassigned cores 2 and 5 are treated as sin- 
gle). This low number of single cores seems highly improbable, 
since there are almost twice as many filaments as cores. A Monte 
Carlo simulation shows indeed that the probability of having 5 
or less single cores in the cloud is less than 1%, and therefore, 



that the observed distribution of cores among filaments does not 
arise from a process equivalent to throwing randomly 19 cores 
into 35 independent filaments. Core formation in L1495/B213, 
therefore, seems to have occurred in a non-random manner, with 
a small fraction of filaments forming multiple cores ("fertile" 
filaments) and the majority of the filaments remaining sterile. 

The dichotomy between "fertile" and "sterile" filaments pro- 
vides a simple explanation for our finding in section |4] that there 
is an excess of cores with separations around 0.25 pc. This is 
so because if most of the cores in the cloud have formed inside 
a small number of fertile filaments, with each fertile filament 
generating more than one core, the resulting core population is 
expected to be distributed in small groups, and each core is ex- 
pected to have a neighbor at a distance on the order of a fraction 
of the filament length. Since the median filament length is ap- 
proximately 0.5 pc, the excess of cores should occur at separa- 
tions around 0.25 pc, which is what is shown in Fig. [7] 

The fertile/sterile dichotomy of filaments, together with the 
multiplicity of cores per fertile filament, points to a core forma- 
tion mechanism that depends more on a global property of the 
filament than on a local event in its interior. To explore this possi- 
bility, we have compared the physical properties of the two fila- 
ment populations using the data in Table [3] One noticeable trend 
is that fertile filaments have a larger mass per unit length than 
sterile filaments: 24 vs 10 M pc" 1 when comparing the means 
and 18 vs 9 M Q pc" 1 when comparing the medians. Although 
the mass per unit length is undoubtedly an uncertain parameter 
due to its dependence on the C 18 abundance, part of this uncer- 
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Fig. 13. Distribution of parameters in the 35 cloud components 
derived using the FIVe algorithm, (a) Aspect ratio, (b) length, 
(c) mass per unit length, (d) non-thermal velocity dispersion, (e) 
global velocity gradients, (f) local velocity gradients. 



tainty is mitigated by the relative nature of the comparison, and 
by the fact that all cores have been analyzed in the same way. 
Interestingly, the data in Table [3] show that fertile filaments not 
only have higher mass per unit length, but that both their mean 
and their median values exceed 16 M pc~', which is the equi- 
librium limit against fragmentation for an isothermal cylinder at 
10 K ( |Stod6Mewicz| 1 963 [ [Ostriker| 1 964} . Sterile filaments, on 
the other hand, have on average a mass per unit length that is 
below the fragmentation limit and therefore should be gravita- 
tionally stable. If this is correct (and it should be tested with a 
more accurate mass estimate, preferable from continuum dust 
emission), it would suggest that core formation in the filaments 
of L1495/B213 simply depends on how much mass the filament 
has been able to accumulate. Most filaments seem to not reach 
the fragmentation limit, and therefore fail to form cores, while a 
selected few do so and produce multiple condensations. 



7.4. Velocity coherence of the filament gas 



In LI 5 17, [Hacar & Tafallaj ( pOTT) found that the internal ve- 
locity field of the cores follows the large-scale velocity gradi- 
ents seen in the filaments, suggesting that core formation has 
occurred inside the filaments with little external interaction (in 




Fig. 14. Velocity profiles of filaments with cores. Blue squares 
are C ls O data and red triangles are N2H + data. 



contradiction with the expectation from some turbulent models, 
see 



Klessen et al. 



2005) . To explore this issue in L1495/B213, 
we present in Fig. 14 velocity profiles for all filaments contain- 
ing dense cores. As can be seen, the gas in the dense cores of 
L1495/B213 (red triangles) presents large-scale velocity fields 
that follow the large-scale gradients of the lower-density mate- 
rial traced in C I8 (blue squares), in good agreement with the 
L1517 results. This suggests that the transition from cloud to 
core conditions in L1495/B213 has also proceeded without a sig- 
nificant change in the gas velocity field. Core formation, there- 
fore, seems to have resulted from a quasi-static process that has 
not decoupled the kinematics of the dense gas from its surround- 
ing material. Internal fragmentation of the filaments, and not 
collisions between gas flows, seems therefore the core-forming 
mechanism operating in L1495/B213. 

Although the velocity field in the cores continues the mo- 
tions of the surrounding material, and this suggests that no kine- 
matic changes occur during core formation, the non-thermal ve- 
locity component of the N2H + lines is on average almost half the 
C 18 value. A similar difference was found in L1517, and could 
in principle result from a lower level of local random motions 
in the dense core material. This would suggest that turbulence 
has been dissipated during core formation. As argued for L1517, 
however, it is also possible that the velocity difference between 
C 18 and N2H + arises from a combination of large-scale gas 
motions and the fact that each molecule traces a different column 
of cloud gas. The N2H + emission is very selective of the dense 
gas due to a combination of excitation and chemistry, and only 
samples a small fraction of the gas along any line of sight (as il- 
lustrated by the sparse distribution of N2H + emission in the maps 
of Fig. [2]). The C 18 emission, on the other hand, is sensitive to 
the extended and lower-density cloud material due to its close- 
to-thermal excitation and high abundance, and it samples a much 
larger column density of gas than N2H + . As seen in Fig. 14 all 



filaments present large scale velocity gradients in the plane of 
the sky, and is therefore likely that similar gradients occur along 
the line of sight. Thus, any C 18 spectrum must contain the con- 
tribution from more parcels of gas moving at different velocities 
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along the line of sight than the equivalent spectrum of N2H + , 
and this will undoubtedly cause a larger non-thermal component. 
This effect is confirmed by the recent numerical simulations of 
core-forming filaments by |Smith et al.| l |2"0T2) . 

The final characteristic of the velocity field that we investi- 
gate is the presence of large-scale organized patterns, often in 
the form of quasi-sinusoidal oscillations (see, e.g., filaments 32 
and 34 in Fig. [I4| . These patterns, also seen in L1517, sug- 
gest that the velocity of the gas in the filaments can be corre- 
lated over scales as large as the filament length, or about 0.5 pc. 
Such a level of organization in the velocity field seems to de- 
viate from the expectation for gas belonging to the turbulent 
cloud regime, which is characterized by a random pattern that 
follows a linewidth-size relation ( |Larson|1981| l. Because of that, 
the L1517 filaments were referred to as "velocity coherent," ex- 
tending the term proposed by Goodman et al. ( 1998) for scales 
of dense cores. As Fig. [14] shows, the same attribute of coher- 
ence seems now applicable to the L1495/B213 filaments, even 
if their non-thermal motions are close to the sound speed. The 
L1495/B213 filaments seem to have also decoupled from the 
general turbulent velocity field of the cloud and become sepa- 
rate entities with a coherent velocity field. 

In L1517, some filaments presented a correlation between 
the oscillations of the velocity field and oscillations in the fil- 
ament density, and from that, it was proposed that the os- 
cillations represent core -forming streaming motions (Hacar & 
|Tafalla||20TT] >. The origin of the oscillations in L1495/B213 is 
less clear due to the more complex velocity patterns and the mix 
of cores at different stages of evolution (i.e, starless and proto- 
stellar) next to each other, which indicates that core-formation in 
the L1495/B213 filaments cannot be treated with the simple per- 
turbation formalism used in the case of L1517. Higher angular 
resolution observations of selected filaments are being planned 
to explore the possible connection between velocity oscillations 
and filament formation. 



7.5. Filament bundles: a scenario of filament-forming 
fragmentation 

The picture that emerges from the previous sections is one of 
hierarchical fragmentation. The L1495/B213 complex (~ 10 pc 
long) seems to have fragmented first into velocity-coherent fila- 
ments of about 0.5 pc in length. Then, some of these filaments 
have further fragmented into cores (< 0.1 pc) in an almost quasi- 
static way. Our discussion so far has concentrated on the sec- 
ond level of fragmentation, the one that breaks up the velocity- 
coherent filaments into star-forming dense cores. This process 
is better constrained by the data because both the "before" (fil- 
ament) and "after" (core) stages can be identified by a different 
molecular tracer. In this section we turn our attention to the first 
level of fragmentation, the one by which the large-scale cloud 
produces the velocity-coherent filaments, which unfortunately is 
less constrained by the data at hand. 

A clue to a possible sequence of filament-forming fragmen- 
tation comes again from the spectra of the B211-B213 region 
shown in Fig. 8] As discussed before, these spectra reveal two 
families of velocity components, one near 5.4 km s _1 and the 
other near 7.0 km s _1 . Each of these families further splits into 
two or three individual components, which according to our 
analysis, represent distinct velocity-coherent filaments. Panel 3 
in Fig. [8] illustrates how the velocity differences between the 
components in each family are smaller than the differences be- 
tween the families. This pattern suggests that there is some type 
of hierarchy in the gas velocity field, in the sense that the com- 
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Fig. 15. Spatial distribution of the nine bundles of filaments iden- 
tified by FIVE. The extent of each bundle is represented by a 
colored background, and the axis of the individual filaments are 
represented by red lines. The grey-scale background is the same 



SPIRE map shown in Fig. 12 



ponents are grouped into families, and that the filaments of a 
given family are more closely connected to each other than to 
filaments of other families. 

The idea that filaments are grouped into families is rein- 
forced by the analysis of the physical and chemical properties 
of the gas. Filaments with velocities close to V L sr = 5.4 km s _1 
(blue-shifted group of components in Fig. [8} are associated with 
intense emission in C 18 and SO, and at the same time have very 
weak N2H + emission. As discussed in Sect. [3] these signatures 
indicate that the gas has an early-type composition, which is ex- 
pected if the material has condensed recently from a more diffuse 
state. The group of components near Vlsr = 7.0 km s , on the 
other hand, is associated with bright N2H + emission and dense 
cores containing YSOs, which are signatures of more chemically 
evolved gas. The 5.4 and 7.0 km s _1 families of filaments, there- 
fore, not only differ in their kinematics, but seem to be at differ- 
ent stages of evolution. This behavior strengthens the idea that 
there is a close connection between the filaments inside a given 
family, while there are noticeable differences between the fami- 
lies themselves. 

To further explore the possibility that the filaments in the 
L1495/B213 complex are organized into families, we have used 
again the FIVe algorithm. This time, we have relaxed its pa- 
rameters so that it can identify structures that are more loosely 
connected than the filaments themselves, but that are still notice- 
ably distinct from each other. Motivated by the differences and 
similarities between the 5.4 and 7.0 km s _1 families, we have set 
the velocity gradient threshold required to define a "friend" to 
5 km s _1 pc~' . Also, to favor structures larger than the filaments, 
we have required a minimum of 30 points to define a distinct 
component and use a box size of 180". The results of this search 
are shown in Fig. 15 where the different families of filaments are 



identified using different colors. As can be seen, the filaments in 
the L1495/B213 complex seem to belong to a small number of 
families, which from now on will be referred to as "bundles" be- 
cause of their thread-like appearance. (The number of bundles 
may be smaller than the 9 shown in Fig. [15] since bundle 3 is 
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likely connected to bundle 6 by low-level emission, [Narayanan 
et al.|2008| ) 

While the search for families using the FIVe algorithm is 
strictly based on the kinematic properties of the gas, the resulting 
bundles shown in Fig.[l5]reproduce remarkably well the division 
of the cloud into regions proposed by Barnard and illustrated in 
the bottom panel of Fig. [2] In Sect. [3] we saw that Barnard's 
sub-division is not simply a morphological one, but that re- 
flects a separation of the gas into regions that have had almost- 
independent contraction histories. The similarity we find now 
between these Barnard regions and the velocity-defined bundles 
suggests that the two structures are either coincident or closely 
related, and therefore are likely to have a common origin. As 
discussed in Sect. [3] the different Barnard regions seem to have 
condensed from the diffuse gas with different time scales, indi- 
cating that a large-scale (pc-sized) fragmentation has occurred 
early on in the history of the cloud. If the velocity-defined bun- 
dles correspond to these fragments, bundle-formation must be 
related to the assemblage of the large-scale filamentary cloud, 
and therefore must have preceded filament formation. 

A piecewise formation picture of L1495/B213 seems at first 
to contradict the high degree of organization suggested by the 
single-filament appearance of the cloud at large (~ 10 pc) scales 
(e.g., Schmalzl et al.|201 0; Palmei rim et al.|2 013). This however 
is not necessarily so if the L1495/B213 complex was formed 
by the convergence of large-scale flows, a mechanism favored 
by a number of theoretical considerations and numerical simu- 
lations (e.g., Elmegreen]|1993 Vazquez-Semadeni et al.| [2006 
Hennebelle et al. 2008 ; Heits ch et al.|20Q8] >. In this converging 



flows scenario, it is unlikely that gas spread over 10 pc linear 
scales can synchronize its convergence without producing pc- 
scale irregularities. For example, perturbations in the gas veloc- 
ity on the order of 10% can easily produce over the expected 
10-20 Myr accumulation time of the cloud (e.g., |Bergin et al. 
2004) the 1-2 Myr differences in the convergence time required 
to explain the L1495/B213 data (Sect. [3]). Thus, although the 
1-2 Myr differences in contraction time have had a notable influ- 
ence in the fragmentation of the cloud into bundles of filaments, 
they still represent a small amount from the point of view of the 
full cloud formation process. 

While the fragmentation of the cloud gas into families of fil- 
aments seems to arise from local differences in the time-scale of 
gas concentration, the internal fragmentation of each family into 
a small group of velocity-coherent filaments seems to be intrin- 
sic to the gas. Even B21 1, the youngest region of the cloud, con- 
sists of multiple velocity-coherent filaments despite not having 
yet formed stars or cores. This early appearance of the filaments 
in the gas suggests that their presence results from physical pro- 
cesses that took place at density regimes lower than those traced 
by our C I8 data (approximately < 10 4 cm -3 ). A number of in- 
stabilities are known to occur at those regimes (e.g., Heitsch et 



al.|2008[ ) and may be responsible for this fragmentation, which 
has a typical scale length of 1 .5-2 pc. Observations of lower den- 
sity gas, like those provided by the more abundant CO isotopes 
( Goldsmith et al.|200 8 ), may hold the key to understand the ori- 
gin of these velocity-coherent filaments. 



8. Summary 

We have studied the 10 pc-long L1495/B213 filamentary com- 
plex in Taurus with the goal of clarifying the process of dense 
core formation. Our data consists of Nyquist-sampled maps in 
lines of C ls O, N2H + , and SO, together with partial mapping in 



the 870 fim and 1200 pm dust continuum. From the analysis of 
these data we have reached the following main conclusions. 
Large-scale properties of L1495/B213. The L1495/B213 com- 
plex appears globally as a single large-scale filament, but differ- 
ences across the complex in the YSO population and chemical 
composition suggest that the gas did not assemble all at once. 
Different regions, approximately coinciding with the condensa- 
tions originally identified by Barnard from optical images, seem 
to have condensed with time scales that differ by up to 1-2 Myr 
(Sect. [3]). Embedded in this molecular complex, there is a popu- 
lation of at least 19 dense cores, some of them starless and some 
of them protostellar. The cores are not distributed uniformly, but 
tend to cluster with favored separations on the order of 0.25 pc 
(Sect. [4}. 

Cloud velocity structure. The C ls O spectra often present mul- 
tiple peaks that indicate the existence of overlapping velocity 
components in the gas (Sect. [5]). We have fitted multiple gaus- 
sians to the C ls O spectra, and found that when the fit param- 
eters are represented in position-position-velocity (PPV) space, 
the velocity components appear as coherent structures that are 
well-separated from each other (Sect. |5j. To disentangle these 
components, we have developed a new algorithm that identifies 
and extracts automatically structures in PPV space. The algo- 
rithm, named FIVe, uses a friends-of-friends approach similar 
to that often employed to identify clusters of galaxies in redshift 
surveys (Sect. |6j. With its help, we have identified 35 separate 
components in L1495/B213. These components are filamentary 
and tend to be aligned with the axis of the large-scale cloud. 
They have typical lengths of 0.5 pc, internal velocity dispersions 
on the order of the sound speed, and coherent velocity fields. 
They also have mass-per-unit-lengths close to the fragmentation 
threshold of an isothermal cylinder at 10 K (Sect. 7.2 1. 
Core formation. Despite the large number of filaments in 
L1495/B213, only a few of them have formed dense cores, and 
the distribution of cores among filaments seems not to be ran- 
dom. A few "fertile" filaments are responsible for most cores 
in the cloud, and this suggests that core formation depends on 
an intrinsic property of each velocity-coherent filament, such as 
the mass per unit length (Sect. 7.3 1. Both fertile and sterile fila- 
ments present coherent, large-scale oscillations in their velocity 
field that suggest that they have decoupled from cloud-wide tur- 
bulent motions (Sect. [7~4| i. When a filament contains a core, the 
continuity between the velocity field of the core gas and that of 
the surrounding less dense material suggests that core formation 
occurs with little effect on the gas kinematics, in contradiction 
with the expectation from models of core formation by gas flow 
collisions (Sect. f7~4] i. 

Hierarchical fragmentation. The 35 filaments of L1495/B213 
seem to be grouped in 9 different families, which we refer to as 
bundles. Filaments inside a given bundle have similar kinematics 
and chemical composition, suggesting that they have a common 
physical origin. This structuring of the gas into bundles, fila- 
ments, and cores suggests that fragmentation in the L1495/B213 
complex has proceeded in a hierarchical manner. First, the cloud 
has fragmented into sub-regions (bundles) that coincide approx- 
imately with the condensations identified originally by Barnard. 
This large-scale fragmentation seems to have resulted from dif- 
ferences in the local time-scale of gas condensation. At later 
times, the sub-regions have fragmented into velocity-coherent 
filaments that have typical sizes of 0.5 pc and velocity disper- 
sions on the order of the sound speed. Finally, a small group 
of fertile velocity-coherent filaments have accumulated enough 
mass to fragment quasi-statically into individual dense cores 
(Sect. [Of. 
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Appendix A: Testing and fine-tuning the Five 
algorithm 

Our analysis of the L1495/B213 complex has used the Friends 
In Velocity (FIVe) algorithm to identify and extract individual 
gas components from the molecular line data. As mentioned in 
Sect. [6] this algorithm is based on the friends-of-friends method, 
which has a well-documented record of reliability in the analysis 
of redshift-survey data over the last three decades (Huc hra &| 
Geller|1982 Berlind et al.|2006) . Its usage with molecular-line 
data, however, exceeds the range of previous applications, and 
for this reason, it is necessary to extend the algorithm testing 
using molecular-cloud conditions. In this appendix, we present 
a series of tests of the FIVe algorithm under conditions similar 
to those of the Taurus cloud to illustrate both the successes and 
limitations of the method. 

The analysis of the FCRAO data from L1495/B213 already 
provides a number of opportunities to test FIVe, and different 
runs of the algorithm were carried out to determine how the 
choice of internal parameters (intensity threshold, maximum ve- 
locity gradient, etc.) affects the output results. It is more in- 
formative, however, to test FIVe using model data that have 
well-defined parameters, since this allows comparing the FlVe- 
derived solution with the original model input. Experience using 
the algorithm shows that most of its behavior can be reproduced 
with a model that consists of two filamentary cloud components. 
This model is simple enough to allow a systematic study of the 
dependence of the results on the choice of internal parameters, 
and thus determine the conditions under which the FIVe algo- 
rithm succeeds in separating cloud components, as well as the 
situations where the algorithm fails by either merging the com- 
ponents, fragmenting them artificially, or simply missing them 
due to low SNR. 

The FIVe algorithm uses both velocity and spatial informa- 
tion to identify cloud components. The role of these two prop- 
erties is almost orthogonal, in the sense that FIVe uses different 
internal parameters to determine the velocity and spatial struc- 
ture of the gas. For this reason, in the following two sections we 
test separately the effect of velocity and space, and we do so us- 
ing two models: (i) two cloud components that overlap spatially 
but differ in velocity, and (ii) two cloud components that have 
the same velocity but differ in position. 

A.1. Overlapping components of different velocity 

The case of two components that overlap in space but have dif- 
ferent central velocities captures a unique feature of FIVe, its 
ability to use velocity information to disentangle cloud structure 
confused in the plane of the sky. Our model for this case as- 
sumes that the cloud is made of two cylindrical gas components 
that resemble the velocity-coherent filaments of L1495/B213. 
Each component has a radial profile of column density parame- 
terized with a gaussian of 150" FWHM, and its length is 840", 
which are typical values of the C 18 emission from L1495/B213 
(Sect. |7.2| >. The emission from these two components has been 
modeled with the CLASS program generating a series of synthetic 
spectra that have the same number of channels and velocity res- 
olution as the FCRAO C 18 O(l-0) data. In each spectrum, the ve- 
locity components are represented by gaussian line profiles that 
have a FWHM of 0.48 km s 1 (average C 18 O(l-0) linewidth) 
and a velocity separation that is kept as a free parameter to ex- 
plore the level of success or failure of the FIVe analysis. To make 
the model more realistic, the intensity of the emission, which 
is 3 K in the filament axis, has been allowed to vary randomly 
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Fig. A.l. Analysis of a model consisting of two components that 
overlap in position but are separated in velocity by 1.25 times 
the line width. Top: Position-position-velocity cube with the re- 
sults of gaussian fits to the spectra. Bottom: Results from the 
FIVe analysis using the favored choice of internal parameters 
(see text). Compare with Fig. 
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by 20% (to simulate the observed non-uniform intensity of the 
data), and random noise has been added to each channel of the 
spectra assuming an rms level of 0.3 K (similar to the noise level 
of our FCRAO observations). 

To simulate the reduction of real data, we have processed the 
model spectra following the same steps as with the observations. 
First, we have fitted the spectra with two gaussian components, 
determining the line center velocity of the components at each 
cloud position. These line center velocities have been used to 
generate PPV diagrams like that shown in Fig. A.l and they 



have been input to the FIVe algorithm to attempt recovering the 
original components. To explore the ability of FIVe recovering 
cloud components with different amounts of velocity overlap, we 
have run a series of models in which the components differ in ve- 
locity by 1.0, 1.25, 1.5 and 2.0 times the spectral line FWHM. 
For each case, we have used different choices of the internal pa- 
rameters of FIVe, such as the number of required neighbors, the 
intensity threshold to select points in the first search for friends, 
and the maximum velocity gradient allowed for defining friends 
(Sect.[6]l. If the FIVe algorithm was able to recover the original 
two components of the model, we have classified the test as a 
success. If not, we have classified the test as a failure. 

presents a summary view of the test results with a 



Fig. A.2 



"traffic light" diagram, which consists of a matrix of pie charts. 
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Fig.A.3. "Traffic-light" diagrams for a model with two cloud components that have the same velocity but are separated in the 
plane of the sky. Left: rates of success (green) and failure (red) for a matrix of separation between components (in units of axial 
FWHM, x-axis) and number of neighbors used in the friend definition (y-axis) Righ: same as in left panel but for models of variable 
intensity threshold (in units of rms). The horizontal boxes enclose the results for the choice of pa rame ters favored by our modeling. 



(All models assume a gradient of 3 km s 
choice of the other panel is used.) 



pc , as recommended by the results shown in Fig. 



A.2 



and for each panel, the best 



Each pie chart shows the fraction of successes (in green) and 
failures (in red) for a series of 10 independent tests of FIVe that 
combine a given value of the velocity separation between the 
filaments (in units of the line FWHM, x-axis) and a value of 
the maximum velocity gradient between the points used in the 
definition of friend (y-axis). Other internal parameters of FIVe, 
like the intensity threshold and the number of points used in the 
definition of friend, have been kept fixed at the values used in 
the analysis of the real data, and which are recommended by the 
models without velocity structure discussed in the next section 
and to which they are more sensitive. 



As can be seen in Fig. A.2 for a velocity separation between 



the components of 1 .0 FWHM, no choice of the velocity gradi- 
ent produces a reasonable (> 50%) fraction of successes. This is 
not surprising, since at this level of proximity, the velocity com- 
ponents in the spectrum start to merge, and the gaussian fitting 
algorithm often converges to a wrong result. At the other extreme 
of separations, velocity components differing by 1.5 FWHM (or 
more) are disentangled successfully by FIVe for any reasonable 
choice of the velocity gradient. This level of success demon- 
strates how FIVe can recognize easily cloud components that 
overlap in space but are separated in velocity by an amount large 
enough to produce a double-peaked profile in the spectrum. 



The critical case in Fig. A.2 is that of velocity components 



separated by 1.25 FWHM (middle column of pie charts). In this 
case, the success rate of FIVe depends slightly on the choice of 
the velocity gradient, and is significantly lower (60%) for the 
smallest value of 1.5 km s _1 pc -1 . Such small velocity gradient 
makes FIVe prone to artificially fragment the components, since 



small errors in the velocity determination (due to the partial ve- 
locity overlap between components) may create gaps between 
sections of the same component that FIVe fails to recognize 
as connected. Although not seen in the simple two-component 
model of Fig. A.2 using a large value for the gradient, like 



5 km s pc , often leads to the opposite effect, the artificial 
merging of the two components that are close in velocity if the 
SNR is low. As a result, the best choice for the maximum veloc- 
ity gradient used in the friend determination seems to be close to 
3.0 km s pc -1 , which is the value we have used in the analysis 
of the real L1495/B213 data. 



A.2. Non overlapping components of the same velocity 

Our second set of models tests the ability of FIVe to separate 
cloud structure based on spatial information. These models as- 
sume that the cloud consists of two components that have the 
same radial velocity but are separated in the plane of the sky by 
an amount that we treat as a free parameter. These models are 
complementary to those of the previous section, and can be used 
to constrain a different set of internal parameters of FIVe, those 
that are related to spatial information, like the number of points 
required to define a friend and the intensity threshold required to 
consider points in the search for friends. 

For this second set of tests, we have followed the same steps 
described in the previous section. First, we have generated a se- 
ries of model spectra using the CLASS program. For this, we have 
assumed that each individual cloud component has the same 
physical dimensions, velocity structure, emerging intensity, and 



20 



A. Hacar et al.: Hierarchical core formation in L1495/B213 




1.25 

5V (FWHM) 



Fig.A.2. "Traffic-light" diagram showing the rates of success 
and failure of the FIVe algorithm when disentangling a model 
cloud that consists of two components that overlap in position 
but differ in velocity. The diagram represents a matrix of pie 
charts, each of them showing the fraction of success (green) and 
failure (red) for 10 independent models with a given velocity 
difference between the components (in units of the spectral line 
FWHM, x-axis) and a maximum velocity gradient used in the 
definition of friend (y-axis). The horizontal box encloses the re- 
sults for the favored gradient of 3 km s~ 1 pc _1 . (All models use a 
number of neighbors of 5 and an intensity threshold of 6 rms, as 
recommended by the results shown in Fig. A. 3 ) 



noise level as the components in the previous section. Then, we 
have fitted gaussians to the spectra, although this time, we have 
used single gaussians because the two cloud components do not 
overlap in space and the spectra consist of single peaks. Finally, 
we have input the fit results into FIVe, and we have compared 
the number and properties of the cloud components found by 
the algorithm with those originally put in the model. As before, 
we have classified the result as a success if FIVe recovers the 
original structure. If not, we have classified it as a failure. 



Fig. A. 3 summarizes the results of our second set of models 
using again "traffic light" diagrams. As mentioned before, the 
models in this section explore mainly two parameters of FIVe, 
the number of points used in the definition of friends and the in- 
tensity threshold used for the selection of points. For this reason, 



Fig. A. 3 consists of two panels, each one showing the success 



rate of one of the parameters when disentangling two compo- 
nents that are separated by a distance value given in terms of the 
filament width (FWHM). 

As can be seen, no choice of parameters allows FIVe to 
disentangle filaments that are closer than 1.25 times the spatial 
width or less. This limiting value is larger than the 1.0 times the 
velocity linewidth found in the previous section, indicating that 
FIVe is not as sensitive when using spatial information as it is 



when using velocity data. The (slight) asymmetry between space 
and velocity likely arises from the simpler treatment that FIVe 
makes of the spatial information, which in turn results from the 
more complex nature of this parameter, that can vary in two di- 
mensions, while the radial velocity has only one dimension to 
change. 



Another result seen in Fig. A. 3 is the success of FIVe to dis- 
entangle components that are separated by 2.0 FWHM (or more) 
in distance unless some of the internal parameters of the algo- 
rithm are set to extreme values (e.g., intensity threshold larger 
than 8 rms). This level of success indicates that FIVe can re- 
construct correctly cloud structure that is well-separated in the 
plane of the sky, and that it only tends to fail when the two cloud 
components overlap in space. 

Fig. A. 3 also indicates that the success rate of FIVe can be 
maximized by choosing appropriate values for both the number 
of neighbors and the intensity threshold used in the friend defi- 
nition. As can be seen, a choice of 5 neighbors and a threshold 
equal to 6 times the intensity rms provides a high success rate 
even in the case of components separated by 1.5 FWHM. For 
this reason, the previous choice of values was used in the anal- 
ysis of the L1495/B213 data. It is important to stress, however, 
that this parameter choice depends on the physical properties 
of the gas and the sampling and sensitivity of the observations. 
Different cloud conditions or a different type of observation may 
require using different parameter choices. 

To finish our analysis, we study how the FIVe algorithm fails 
in the limiting case of components separated by 1.5 FWHM 
when its parameters have not been optimized. This is instruc- 
tive because it shows the compromise that is often needed to 
determine the optimal value of a given parameter. As Fig. A. 3 



shows, in the case of a variable number of neighbors (left panel), 
a choice of 5 neighbors produces a higher rate of success than a 
choice of 3 or 7 neighbors. This is because a choice of a too 
few neighbors, like 3, can lead to the erroneous merging of the 
components, since it becomes possible that a statistical fluctua- 
tion due to noise or variations in intensity affects a few positions 
that can find enough neighbors to create an artificial bridge be- 
tween separated components. On the other hand, a choice of too 
many neighbors, like 7, can lead to the erroneous splitting of a 
single component. This is because a requirement of 7 neighbors 
for the definition of a friend leaves too little room for possible 
noise effects. A statistical fluctuation that eliminates several of 
the 8 possible nearest neighbors breaks the component in two 
and makes the algorithm fail. 

A similar situation occurs with the choice of the intensity 
threshold, as shown in the right panel of Fig. A. 3 Choosing too 
low a value (like 4 rms) allows the low-intensity tails of the gaus- 
sian distributions to merge, and this makes the algorithm fail to 
separate the two cloud components. Choosing too large a thresh- 
old (like 8 rms), on the other hand, makes the algorithm sensitive 
to the loss of points in statistical fluctuations, and this leads again 
to artificial fragmentation. The optimal value of choice (thresh- 
old of 6 rms), therefore, represents a compromise between the 
two opposite trends of merging and fragmentation, and makes 
it possible to disentangle cloud components separated by 1.5 
FWHM. 

A.3. Summary: successes and limitations of FIVe 

While our two-component cloud can only be considered a toy 
model, its analysis with FIVe presents remarkable similarities 
with the analysis of the real L1495/B213 data. This is likely a 
combination of the realistic emission parameters of the model 
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Fig. B.l. Schematic view of the two-step process used to deter- 
mine the principal axis a filamentary velocity component. The 
grey-scale and contours represent the C 18 O(l-0) emission. The 
blue straight line is the first-step determination of the principal 
axis, and the red polygonal curve is the second-step determina- 
tion of the axis. 



and the fact that most problems analyzing of real data arise when 
dealing with pairs of components, especially when trying to sep- 
arate them if they are close in velocity or position. Under these 
conditions, the previous tests show that FIVe can recognize com- 
ponents with an accuracy that approximates that of a human ob- 
server who has the patience of inspecting the individual spec- 
tra and who looks for connections between spectral features at 
different positions. Reproducing such a human behavior was in 
fact the main goal when designing FIVe, given the large size of 
the L1495/B213 dataset and the need to automatize the velocity 
analysis. 

The success of FIVe is remarkable in view of its simple- 
minded analysis, but is accompanied by important intrinsic lim- 
itations. As we have seen, FIVe cannot disentangle components 
that are blended either in position or velocity by amounts com- 
parable to their spatial or velocity width, even if the blending oc- 
curs over a limited region in space. Improving on that limit may 
require a significantly different type of analysis. Like similar al- 
gorithms, FIVe has only a local view of the emission, and cannot 
take advantage of the information provided by the data on global 
scales to separate so-far unresolvable components. Adding such 
a global view to FIVe not only requires a more sophisticated 
numerical scheme to combine local and global patterns, but an 
a-priori understanding of the physical nature of the components 
themselves. This understanding, unfortunately, still eludes us to- 
day. 

Appendix B: Principal axis determination of the 
velocity components 

To analyze different properties of the filamentary velocity com- 
ponents identified with the FIVe algorithm, such as their spatial 
location and velocity field, it is convenient to determine for each 
of them an axis that follows the long dimension of the filament. 
We have done this using the process illustrated in Fig. |B.l| which 
consists of two steps. First, we have fitted a straight line to the 
distribution of all points identified by FIVe as belonging to the 
same velocity component, weighting each point by the intensity 



of the C 18 O(l-0) emission (in cases of suspected CO depletion, 
the C ls O intensity has been corrected as described in Sect. ^J, 
As a result of this step, we have obtained a first approximation 
to the principal axis of each velocity component (blue line in 
Fig. - 



B.l 



Using this axis as a reference, we have then divided 
the set of points into 90"-long segments, and for each segment, 
we have determined the emission centroid following the same 
wighting scheme as in the first step. Connecting the emission 
centroids with a polygonal curve, we have defined the second 
and final approximation to the principal axis of each filament 



(red line in Fig. B.l 
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